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Abstract. Problem of meteor orbit determination for a given parent body is discussed. 
Some of the published methods for obtaining meteoroid's orbital elements at the moment 
of intersecting Earth's orbit on the basis of geometrical variation of parent body's orbital 
elements are discussed. The main result concerns the following two facts: i) in real situ- 
ations physical quantities for the change of orbital elements of the parent body must be 
used, and, ii) the usage of Southworth and Hawkins (1963) D-criterion yields results not 
corresponding to observations. 
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1. Introduction 

In searching for the relation between parent body (comet, asteroid) and its possible 
meteor stream, one needs to calculate theoretical orbit of a meteoroid at the time when 
it may become a meteor. Several methods have been worked out during the last several 
tens of years. Their summary may be found in Svoren et al. (1993). Our aim is to discuss 
some of the published methods and improve them. We also suggest physical method. 

Svoren et al. (1993), Neslusan et al. (1994), as examples, use D-criterion of Southworth 
and Hawkins (1963) in many of the methods. Moreover, in final comparison between the- 
oretical and observational result, this D-criterion is considered to be a decisive quantity. 
However, it seems to be not adequate to use this criterion due to its many "curious" 
properties (Klacka 1995). 



We will use standard orbital elements in the paper: q perihelion distance, e - eccen- 
tricity, i - inclination, Q longitude of the ascending node, lo - argument of perihelion. 
Moreover, subscript "P" refers to parent body, subscript "M" to meteoroid which inter- 
sects the orbit of the Earth. The Earth's orbit is supposed to be circular. 

2. General considerations 

If the meteoroid intersects the Earth's orbit, then one node of the meteoroid's orbit is at 
distance r^ — 1 AU from the Sun. Thus, we have (/ - true anomaly) 

sin {ljm + ./m) = and 

TE = — -r- ■ (1) 

I + CM cos/m 

Thus, if the ascending node (s > => sin(ajM + /m + e) > - this corresponds to 
^M + Jm ~ 2 tt) intersects the orbit of the Earth, then 

qm (1 + Cm) 



lAU = rAE = 



(2) 



1 + CM COSO^M 

If the descending node (e > =^ sin(cjM + /m + e) < - this corresponds to 
^M + ./m = tt) intersects the orbit of the Earth, then 
qm (1 + Cm) 



1 AU = roE 



1 



CM COSWM 



(3) 



3. Methods of variation of the orbit in the orbital plane 

Three orbital elements may vary: q, e, lu. If lum = ^p and qM = qp, or, lum = tjp and 
sm — ep, the situation is trivial - in both cases only one orbital element (e, resp. q) is 
varied in order to Eq. (2) or Eq. (3) be fulfilled. 



3.1. The argument of perihelion is varied 

Let us suppose that it is sufficient to make a change of LUp, only, for the purpose of that 
a meteoroid may intersect the Earth's orbit: qp = qM, ep = cm- Eq. (2) yields then 



(4) 



wmi — arccos < 


r 1 

[ep 


■ gp (1 + ep) ^ 

TAE 


WM2 = 360° - 


^Ml • 


Eq. (3) yields 




wm3 — arccos < 


r 1 

[ep 


' ^ gp(l + ep)' 
roE 


<.^M4 = 360° — 


WM3 





(5) 

If ll!m does not exist, then the assumption that the change of only one element ujp is 
sufficient in obtaining meteor, is incorrect. 



3.2. Perihelion distance and eccentricity are simultaneously varied 

Let us consider now that lum = ^p s^nd that q and e may change. If we use also D- 
critcrion of Southworth and Hawkins (1963), as it is done in Svorch et al. (1993) (only q 
and e are different for parent body and meteoroid) , the minimum simultaneously change 
of q and e is given by the condition that the function 

^Ih = ilM - qpf + (eM - epf (6) 

gains its minimum value. Besides this condition, Eqs. (2)-(3) also hold [ujm = ^p)- If 
we substitute qM from Eq. (2) into Eq. (6), we obtain Dsh as a function of cm only and 
its extremal - minimum value can be easily found. The same procedure must be replied 
for Eqs. (3) and (6). In both cases it can be analytically proved that only one solution 
corresponds to cm > 0. For some parent bodies the minimum value of cm corresponds 
to CM G (0, 1), for the others cm <= [1,9/4) (the last value is approximate). In any case, 
the value of Dsh obtains the value, which is even in the interval cm G (0, 1] less than it 
would be if only one of the elements qM and cm is varied. 

These analytical results are in contradiction with the results presented in Svoren et 
al. (1993 - see Tables 1, 3, 4, 6). 

Another important result shows Table 5 in Svoren et al. (1993): the simultaneous 
optimal change of qM and cm yields value of Dsh which is less than the value of Dsh 
if only qM is changing, but the concidence between theoretical calculations and observed 
data is better in the latter case. This result confirms the fact that D-criterion of South- 
worth and Hawkins is not good approximation to the real processes, as it is discussed in 
Klacka (1995) from the general point of view. 

4. Methods of variation and rotation of the orbit in space 

We will discuss only one method here, in order to stress again that we do not need any 
type of D-criterion. 

^.1. Adjustment of the orbit by rotation around the line of apsides 

We will suppose qM = qp, cm ~ ep. 

Let the coordinate system S' be created in the form that the orbital plane of the parent 
body is characterized by the condition z' = and let the perihelion lies on the positive 
part of the .t'— axis: perihelion is characterized by the unit vector x' = 1, y' = z' = 0. 
This unit vector has coordinates x, y, z in the original coordinate system S (ecliptical 



system, plane of the ecliptic: z — 0): 

Xpp ~ cosflp cos Lup — sin rip sincjp coszp 
VpP = sinilp coso^p + cosilp sincjp cosip 
Zpp = sinwp sinzp . (7) 

Let the unit vector normal to the orbital plane of the parent body has coordinates x^ = 
y^ = 0, z^ = 1. This normal unit vector is characterized by the following coordinates in 
the system S: 

XnP = sinilp sinip 

ynp = — cos Op sin ip 

ZnP — cosip . (8) 

Now, we make a rotation around the x'— axis in an angle <I>. As a result we obtain a 
new orbit - orbital plane of the meteoroid - characterized by new normal unit vector, 
which has coordinates in the system S: 

XnM = sin Om sin Im = 

= (cos Op sincjp + sin Op coso^p cosip) sin$ + 
+ (sin Op sinzp) cos$ 

VnM = - cosOm sinzM = 

= (sin Op sino^p — cos Op coso^p cosip) sin$ — 
— (cos Op sinip) cos$ 

ZnM = COSiAf = 

= — (coswp sinip) sin(f> + (cosip) cos$ . (9) 

Unit vector of the perihelion of the meteoroid's orbit is 

XpM — cosO^f COS lum — sIuOm sint^Af cos iM 
VpM = sIuOm coswm + cosOm sino^M cos iM 
ZpM ~ sincjM sin iM , (10) 

and it is the same unit vector as the unit perihelion vector of the parent body (x'— axis 
- the line of apsides - is fixed during the rotation) . 

Thus, we have one set of three equations given by the equality of the right-hand-sides 
of Eqs. (7) and (10) and another set of three equations given by Eq. (9). For a given angle 



$ Eq. (9) determines quantities im and JIm, and, then, Eqs. (7) and (10) determine the 
angle ojm- The angle $, for which Eqs. (2) or (3) hold, may be easily found. 

We have found the angle of rotation $ around the line of apsids which determines 
the orbit of a meteoroid to become a crosser of the Earth's orbit. Thus, we have found 
all the required elements: im, ^m and ujm- We do not need any knowledge about the 
D-criterion of Southworth and Hawkins as it is in the method presented in Svoreh et al. 
(1993). 

5. Physical Methods 

Let us consider physical access to the problem. 

5.1. Energy-Hz-TT Method 

By this method we consider simple model in which the following quantities of the me- 
teoroids are equivalent to those of the parent body: semimajor-axis a, z-component of 



angular momentum -^/a (1 — e^) cosi, tt ^ lu + fl. Considering Eqs. (2) and (3), the 
minimization of the magnitude of the difference of the angular momenta vectors between 
meteoroid and its parent body |Aif | = \H — Hc\ (global and local minima) yields 
several theoretical radiants - those with the smallest values of | AiJ| should be realized. 

5.2. Energy-Ci-C2-C3 Method 

By this method we consider the model of Babadzhanov and Obrubov (1987) in 
which the four quantities of the meteoroids are equivalent to those of parent body: 



semimajor-axis a, z-component of angular momentum y^a (1 — e^) cosi = C\, C2 — 
e^ (0.4 — (sinz)^ (sino;)^), C3 = tt = lu + il. Possible orbits of meteoroid and the cor- 
responding radiants are found on the basis of these four conservation quantities. The 
physical idea of the method yields that all the obtained radiants are equivalently signifi- 
cant and they all should exist in reality. 

6. AppHcation of the Physical Methods and Method of the Section 4.1 

As an application of the mentioned physical methods we present theoretical orbital ele- 
ments and radiants of meteoroids of the comet 73P/Schwassmann-Wachmann 3 for the 
year 2006. Orbital elements of the comet, considered in our calculations, are: Qc ~ 0.93973 
AU, ec = 0.69331, ic = 11.398°, lUc = 198.824°, n^ = 69.883°; these elements were ob- 
tained by the authors of the paper Gajdos et al. (1998), where also theoretical radiants 
by other methods (sections 3 and 4 of this paper) are presented. 



The Encrgy-H^-TT method yields, in principle, four theoretical radiants. The values of 
the last column of the Table 1 suggest, however, that only one radiant should correspond 
to reality - that given by the first row in Table 1. 

Table 1. Energy-Hz-7r Method. 

q\Av\ e in un nn x^n an sn lAgp 

0.9470 0.6910 12.3 209.5 59.2 59.2 211.4 26.9 0.0009 

0.9176 0.7001 8.1 143.1 125.6 125.6 172.0 28.9 0.0143 

0.8954 0.7078 0.1 318.4 310.3 130.3 163.5 6.7 0.0205 

0.8954 0.7078 0.1 41.6 227.1 47.1 193.9 - 6.2 0.0207 



The Energy-Ci-C2-C3 method yields also four radiants. The physical idea of the 
method predicts the real existence of all the four radiants (see Table 2). (We mention 
that no meteor stream existed for the year 1930 according to Energy-Ci-C2-C3 method.) 



Table 2. Energy-Ci-Ca-Cs Method. 



q[AU] e z[°] c.[°] »[°] AqP] «["] 8n 

0.9263 0.6977 9.6 214.9 53.9 53.9 206.7 18.8 

0.9263 0.6977 9.6 145.2 123.6 123.6 173.9 33.1 

0.9263 0.6977 9.6 325.2 303.6 123.6 151.1 - 17.9 

0.9263 0.6977 9.6 34.9 233.9 53.9 183.7 - 32.1 



Finally, we present also radiants for the year 2006 obtained by the method described 
in section 4.1. Again, also this method yields, in principle, four theoretical radiants. The 
values of the last column of the Table 3 suggest, however, that only one radiant should 
correspond to reality - that given by the first row in Table 3. (We mention that no meteor 
stream existed for the year 1930 according to this method - Qc > 1.0 AU.) 

The first two rows of Table 2 correspond to the first two rows of Table 1 . The third 
row of Table 2 corresponds to the second row of Table 3. However, the last columns of 
Tables 1 and 3 show that only the first row of Tables 1 and 3 can correspond to the reality. 
In other words, the real existence of all the four radiants given by the Energy-Ci-C2-C3 
method seems to be improbable. 



Table 3. Adjustment of the orbit by rotation around the hue of apsides method - method 
decribcd in section 4.1 

4AU] e in ojn »[°] AeP] aU Jp] <i>[°] 

0.9397 0.6933 7.0 211.5 57.1 57.1 203.1 14.0 4.8 

0.9397 0.6933 7.0 328.5 299.7 119.7 154.0 - 12.9 16.8 

0.9397 0.6933 173.0 328.5 237.1 57.1 337.7 - 13.5 184.8 

0.9397 0.6933 173.0 211.5 119.7 119.7 19.2 12.3 196.8 

7. Discussion 

We have discussed some of the pubUshcd methods of calculating meteoroid orbits crossing 
the Earth's orbit if the parent body's orbit is known. 

If the parent body's orbit is changed in so a trivial geometric way as it is often done, 
one should try to avoid to any further requirement, e. g., to use D-criterion of Southworth 
and Hawkins. The reason is that this criterion yields nonphysical results, as it is pointed 
out in the last paragraph of the section 3.2 and discussed in Klacka (1995). 

The use of simple methods has one great advantage - calculations are very simple. 
However, since changes of orbital elements in such simple methods are only of geometrical 
character, the methods do not correspond to any real physical process - the methods do 
not use any physical basis. Since real physics is very complicated and we do not know 
exact perturbations on mcteoroids (the importance of nongravitational effects), we have 
to make predictions for a new possible pairs "parent body - meteor stream" on the basis 
of a simple physics, or, on the basis of the known pairs. 

Possible way is to find a function F (containing five orbital elements in an independent 
way) which is a good approximation to reality and does not exhibit any inconsistencies. 
Suggestion of this type may be found in Klacka (1995), where also general method is 
presented - special case of this method is Energy-Hx-Tr method. The method in reality 
corresponds to the following set of equations: 

Am = fi{(3jp) , i,j^lto5, (11) 

where /?— s represent orbital elements (P - parent body, M - meteoroid, meteor 
stream) and Eqs. (2) or (3) hold. The unknown functions must be found (approximated) 
on the basis of observational data - pairs "parent body - meteor stream" . The conditions 
fi = PiP {i = 1 to 5) hold if /3jp fulfill Eqs. (2) or (3). Thus, the method presented in 



Klacka (1995) is equivalent to the following sets of equations: 

Am = P^p + l^^iijL^^L _ A f^^^p ) , or 
11 + ep cos ojp J 

Am = Ap + |^dl_t^£L _ A /2^(/j ) , 
[^ 1 — ep cos bjp J 

i,j = lto5, (12) 

where fu and /2i are finite, and, moreover, Eqs. (2) or (3) hold. 

In any case, one must be aware of the fact that real meteoroid streams have various 
dispersions of orbital elements and that our methods of calculating new pairs "parent 
body - meteor stream" are only a first approximation. 

In principle, no reasonable metric for various trajectories of parent bodies 
and their meteoroids exists. Gravitational and nongravitational forces acting 
on parent bodies and their meteoroids are different in various parts of orbital 
element's phase space, and, no global metric can be constructed on the basis of 
the knov^r pairs "parent body meteor stream". Useful methods for predicting 
new pairs "parent body meteor stream" must be based on the known pairs 
"parent body meteor stream" and inevitable weighting factors, functions of 
parent body's orbital elements, present in "metric" disturb the basic property 
of metric its symmetricity. 

It is hard to make reasonable predictions for meteor streams corresponding to new 
parent bodies lying outside the zone of the known parent bodies for which we know the 
relation "parent body - meteor stream" . Only extrapolations can be made in this case. 
Much better predictions are obtained for the inner part of the zone of the phase space. 
Energy-H^-TT method seems to yield good results in any case. 

8. Appendix 

Derivation of equations of section 4.1 is presented here, on the request of the referee. 

Let the Cartesian coordinate system S' is created from the Cartesian coordinate sys- 
tem S by the rotation characterized with Eulerian angles fi, i and lu. If x, y, z are 
coordinates of a vector in the system S, then x' , y', z' are coordinates of the same vector 
in the system S'. We have 

fx\ /Cii Ci2 Ci3\ fx'\ 

y = C21 C22 C23 y' 

\z) VC3I C'32 C33; \z' j 

and the transformation matrix C is of the form 



COS ft COS uj — sin 51 sin oj cos i — cos 17 sin ut — sin 57 cos oj cos i sin 57 sin i \ 

sin 57 cos uj + cos 57 sin uj cos i — sin 57 sin uj + cos 57 cos uj cos i — cos 57 sin i 

\ sin a; sini cos a; sini cosi / 

(Transformation is orthogonal and so inverse transformation can be easily obtained, using 

transposed matrix.) 

Unit vectors of perihelia are defined by conditions x' = 1, y' = z' = 0, and so the last 

equation immediately yields Eqs. (7) and (10). Unit vector normal to the orbital plane 

is defined by conditions x' = y' = 0, z' = 1, and, again, the last equation immediately 

yields Eq. (8). 

Moreover, if we make a rotation characterized by an angle <& around a;'— axis, we 

obtain 



y" = y' cos$ + z' sin$ 

z" = - y' sin $ + z' cos $ . 

Now, if we put x" = y" = 0, z" = 1, then x' = 0, y' = — sin$, z' = cos$, and the first 

equation of the appendix yields results which are consistent with Eq. (9). 
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